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Abstract. Two- and three- particles photodisintegration of the triton is investigated in a three-dimensional 
(3D) Faddeev approach. For this purpose the Jacobi momentum vectors for three particles system and spin- 
isospin quantum numbers of the individual nucleons are considered. Based on this picture the three-nucleon 
Faddeev integral equations with the two-nucleon interaction are formulated without employing the partial 
wave decomposition. The single nucleon current as well as tt— and p— like exchange currents are used in 
an appropriate form to be employed in 3D approach. The exchange currents are derived from AV18 NN 
force. The two-body t-matrix, Deuteron and Triton wave functions are calculated in the 3D approach by 
using AV18 potential. Benchmarks are presented to compare the total cross section for the two- and three- 
particles photodisintegration in the range of < 30 MeV . The 3D Faddeev approach shows promising 
results. 
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1 Introduction 



and the necessity of complicated recoupling algebra that 
accompanies partial wave based calculations. 



The study of 3N system with the electromagnetic inter- 
action shed light on the ambiguity of different kinds of 
nuclear forces participating in the few body systems. Fol- 
lowing the introduction of Faddeev formulation for the 
three-body system [l]-[2], a new effort using this scheme 
was started. The electrodisintegration [B] and photodisin- 
tegration [3] of '^le and ^ in the Faddeev scheme were 
performed as early attempts in this respect. The photo- 
disintegration calculation based on Faddeev formalism for 
3N bound and continuum with the same 3N Hamiltonian 
were performed with quite simple separable NN interac- 
tions [S] . The other approaches in this field can be consid- 
ered to be Green Function Monte-Carlo f6^ and Lorentz 
Integral Transform [7] methods. 

The photodisintegration of "^le was investigated by 
pair-correlated hyperspherical harmonics method using 
AV18 NN and Urbana 3N forces [5]. Also, a systematic 
application of Faddeev formalism was introduced, based 
on the partial wave decomposition, for the calculation of 
the photodisintegration of ^le and '^I with the two- and 
three- body forces [S]- 

We use the novel three-dimensional (3D) approach pHl 
Pl|[ni|T?lPllP^[rei[r71[TFl[TOl[^ ■ In this approach we em- 
ploy the momentum-vector variables as basis states. We 
derive the Faddeev integral equation in a realistic 3D scheme 
as a function of Jacobi momentum vectors and the spin- 
isospin quantum numbers. We implement this 3D scheme 
in the two- (Nd) and three-body (3N) photodisintegra- 
tion of ^H. This 3D approach avoids truncation problems 



Nd 



It is our aim to calculate the nuclear matrix element 
.ch is given by A^^Af ^ i i^) {(j)„\p\u) or 

_ for j^H -> 3A^ or -/^H Nd, re- 

spectively. We introduce '^-'(0o| and {(f>d\ as outgoing 3N 
and Nd system, P as permutation operator and \U) — 
{1 + P)0\>1') + {tGoP . . . ) I [/) as an auxiliary state with O 
as the current operator. The wave function \l/ is the three- 
body bound state, t is the two body transition operator 
and Go is free propagator. The second part oi\U) is related 
to the solution of | U) of this equation. We introduce a for- 
malism to solve this equation directly in the momentum- 
vector variables avoiding the partial wave (PW) decom- 
position scheme. 

We use the AV18 potential for the calculation of wave 
functions of the triton and deuteron as well as the two- 
body t-matrix. For nuclear electromagnetic current oper- 
ator we choose, in addition to the single-nucleon current, 
the two-body contribution in the form of the tt— and p— 
meson exchanges with connection to NN force AV18. 

This manuscript is organized as follow: In section[5] for- 
mulation of nuclear matrix equation in the 3D approach 
is introduced and |C/), as an auxiliary state, is derived 
for the two-body interaction with handling of singular- 
ity problem. In this section the appropriate formulation 
is constructed for the one- and two-body currents to be 
implemented in the 3D formalism. In section [3] numerical 
method and results are presented and section 2] is con- 
cluded by a summary and outlook. 
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2 Formulation of N matrix equation in a 3D 
approach 

Our aim is to calculate the total cross section of the triton 
photodisintegration. Final states can be either Nd or free 
3N state. For each we can write the cross section as a 
function of the tensor components of the nuclear matrix 
element, which is the nuclear current sandwiched between 
initial bound and final scattering states: 



Nd 



2Q \±n_ _ kd k,. 



4 E (Ki 



Nd\2 



\N 



Nd\2\ 



(1) 



SAT 



1 



dkidk2dki^ ^ -— ( 

m,mi,m2:m3 



kik2mN 



I Q cos 6*2 — ki cos 6*12 — 2^2 



(2) 



In the above equations tojv is the nucleon mass, a w 
is the fine-structure constant, k„ and k^ are the neu- 
tron and deuteron momentums, respectively, k^s are free 
particle momenta and 9iS are their angles with Q as the 
momentum of the photon. O12 is the angle between mo- 
menta fci and /c2- There is also a summation over all the 
spins of ingoing and outgoing particles which are indicated 
by TO, TOd, TOi, m2 and TO3. 

To obtain the nuclear matrix element we follow the 
formalism in 



N 



3N 



1 



(3) 



N"-^-{MP\U), (4) 
where |(/>o)^~'' is the three particles scattering state as; 

|0o)<-^ = (l + Goi)|</)o), (5) 

so we can write: 

N'"" = ^[{cl>o\P\U) + {cl,o\tG,P\U)] (6) 
and is Nd state 

= |qTOii/i)|V'd), (7) 

where tpd is the deuteron wave function, and \U) is the 
auxiliary state that can be calculated by solving the fol- 
lowing equation. 



\U) - {l + P)J\ip) +tGoP\U). 



(8) 



1-0) is the three-body bound state, J is the current. Go 
is the 3N free propagator, t is the two-body transition 
operator and P is the permutation operator. 



2.1 Introducing free basis states 

We introduce |0o) as our free basis states where the par- 
ticles 2 and 3 are in subsystem and particle 1 is the spec- 
tator. The state is antisymmetric under permutation of 
subsystem particles. 



|(/)o) = |pm2TO3i/2i^3)°|q'71ll^l) 



(9) 



where rrii and are spin and isospin of the individual 
nucleons and p and q are Jacobi momentum for a three- 
particle system. 



P= i(k2-k3), 
q=|[ki-i(k2 + k3)]. 



(10) 
(11) 



where ki, k2, k3 are momentums of the particles. So we 
have: 

I pqm 1 m2 TO3 :/i :/2 ^^3 ) ° 
= ^(|pqTOim2TO3i^ii/2J^3) 

-| - pqTOlTO3TO2J^lJ^3i^2))- (12) 



2.2 Nuclear matrix in free basis states 

We now apply these free basis states to the left hand side 
of each term of the Eq. (|S]): 

((/)o|C/) = "(pqTOiTO2m3i/iz/2i^3|C^) 
= "(pqmiTO2m3Z/i;/2i^3|(l + 

-h''(pqTOlTO2TO3Z/i:/2J^3|tGoP|t/). (13) 

In the next step we need to know the effect of the permu- 
tation operator on the free basis states. 

P|pqTOim2TO3Z^ii/2J^3)° 
13 1 
= \i-2P ~ i^)(P ~ 2'^'m-2.rrizmxV2V-iVxT 

13 1 
+ [(-2? + 4q)("P " -j^^in-imxm2V-iVxV2Y , (14) 

so the first part of Eq. © is 

(0O|P|C/) = "(pqTOlTO2TO3i^ll^2J^3|P|C/) 

13 1 
= ''((-2?- ^q)>(p- -j<^m2m-imxV2V-iVY\V) 

13 1 
+ "((-2P+ 4^^'*^^P^ -q)m3mim2i/3i'ii^2|C/)-(15) 

In the evaluation of the first term of Ec^. (jl3p we can write 

°(pqmiTO2TO3Z^ll'2i^3|(l + P)J\^) 

= i^pi\mxm2m-iVxV2V-i\J\^) 

13 1 
+ ''((-2P~ 4q)(P~ 2^)™2"^3"^l^^2^^3^^lkl'^) 

13 1 
+ "((-2?+ ^q)(-p- -q)TO3miTO2i^3;^ii^2k|f)-(16) 
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Evaluation of each term of Eq. p6p is similar and repre- 
sents generally, by using completeness relation and insert- 
ing the free basis states in a suitable position: 

"■{pqmim2m3iyiiy2i'3\J\W) = ^ J dp'dq' 



X ° (pqm 1 TO2 "^3 i^i 1^3 1 •/ 1 p'q' m2 f 2 ) ° 



(17) 



triton wave function in the free basis states is introduced 
symbolically as follows: 

'^m\"^?™3(P'q)= ^V^mim^m^v^v^i^^]^). (18) 
For the second part of Eq. we have: 



E 



dp'rfq'dp"rfq" 



° (pqmi 77127713 vi V2 1'a I i Go | p'q'm'^ mjmg v[ t^2 ^3 ) ° 
x"(p'q'm;m>^i.;i.>^|P|p"q"m'/m>J,V;V^'^r 

a l„" „ll ^11 ^11 ^11 , ,11 , ,11 , ,11 



X»(p"q"m'/m>^V;V^'^'|C/), 



(19) 



there are two important terms in Eq. (IT^ to be evaluated. 
The first term including tGo operators can be written in 
terms of two body t-matrices in the free basis-states mul- 
tiplied by free propagator, Go = "e^~h~i in the energy of 
subsystem. For the second term of Eq. (IT9|) . one can pro- 
ceed to the final form by considering the properties of the 
permutation operator, P ( Eq. (I14p ). and the symmetry 
properties of the free basis states. The final form of the 
second part of Eq. is evaluated as follows: 

{(j}o\tGop\U) = "(pqmi?772m3i/iJ/2i'3|iGoP|C/) 

1 



2 



4m 



-S{q' ~ q) 



X °(pm2m3i/2J^3|i| — n2m27n3t'2t^3)'' 

X (niq"7)73TOl7)72t'3'^l i^2 I ^) I 



(20) 



where the shifted argument IIi and 112 are defined as 
follows 



Hi 



1 



1 



2" ' " ' 2" 

We write for Eq. ^ with the help of Eq. (PH]): 

" (pqm 1 7772 7773 1'l ^2 1's I C^) 

=" (pq777 l7n2r7l3l/ll/2l'3|(l + P)J\il') 

1 



(21) 



+2 E / 

mLmLyL v'. 



dq'' 



X " {pm2m-iV2Vj,\t\ - q' 
1 



^ _ g^+g"^+q. q" 

„ 1 



-q, 777^77737^2^^3)° 



q" , q"7(7^j777 l 777'2 T^gT^l 7^^ | C/) , 



(22) 



where the first term has been described before in Eq. (jl6p . 

Finally the calculation of N'^^ after inserting Eqs. ([T5|) 
and (1201) in Eq. (g]) leads to: 

13 1 



E 



dq"- 



1 



g2_|_q»2^.q. q» 



E 



X ''(P7772 7773 7/2 1^3 1^1 - Q" - ^Wi'^WiT 

X " (q + ^q" , q"™ Wi "^2 ^'3i^i;^2 1 , 



(23) 



we proceed to N^"^ calculation by inserting the complete- 
ness relation of the free basis states and Ec^. dTSl) into: 



X " {p'c^m[ m'^ m'^ v[v'^v'^\P\U) 

= \ j dp'{^,^'^\p'm'2m'A^T 

13 1 

><["(- 2?' - 4^'? - 2^' ^'2'm'z'mivW'i^i\U) 

13 1 

+"(-2?' + 4^' - 2^' "^3'^l"^2^^^l'^2|C^>]- 



(24) 



2.3 Singularity problem and rewriting Eq. (|22[ | 

We know that the two-body t-matrix has a simple singu- 
larity at = -Erf, which is the deuteron binding energy. 
There is also a moving singularity in Eq. (I20p . which is 
very difficult to handle. To solve this moving singularity 
problem we use the method introduced in Ref. [H] and we 
separate the radial part from the angle part in the Dirac 
delta functions. 

,5(p'-i-n2)^(p"-ni) 

then we write: 

° (pqrril 7712 7773 7/1 1/2 7/3 1 tGo P 1 1/) 

=2 E jdc('dp'dp'T-:zr^^,{v>.p'n2) 



xS{p' - n2)S{p" - n^) 



4m m 



(26) 



in the above equation we used: 



^a,.2.3.2.3 ( ) 



" (pr772 7r73 7/2 7/3 I i |p'7T2 7772 7773 7/2 7/3) " , 



(27) 
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and 

= ''(p{-ni)qmim2m3Viiy2i^3\U), (28) 

we have for the radial delta function: 

6{p' - n2)S{p" - TTi) 
1 



(29) 



(5(p'-^/-q2+g"2 + gg"a;") 



we can write by using the property of the delta function: 



6{p' - + q.2 + qq.^. ) ^ ^^(^" _ ^p)^ (30) 



where 



Xo 



qq" 



(31) 



Xo is cos 9qii and it is expressed in this interval: —l<xo< 
+1. 

We use the Eqs. ^ and ^ in Eq. therefore, 
the final form of Eq. (|22|) can be rewritten as follows: 



" (pqmi m2TO3 1'l ^^2 1 1^) 
="■ (pqmim2m3i^iz^2i^3|(l + P)J\'ip) 
4 1 ^ f°° . . p' 



E / dp 



1 /.27r 



E-^-I^ 

4m m 



rfq"?" / dxq>> / d(/)5"5(a;"-2;o) 
'If-p'l J-i Jo 

^^a,v,.W.,f, (p p'Tj )[/-3-i-2 (p"(-n^)^ q"), (32) 



where 



3q"2 3q2 



4 ' 



(33) 



we have used the following relation to extract the singu- 
larity from two-body t-matrix: 



m = 



(34) 



2.4 Current 

The current is consisted of single-nucleon and two-nucleon 
currents 

J =J' + J", 
J' = j\l) + j\2) + j\-i), 
J" = J" (1,2) + J" (2,-3) + J" (3,1). (35) 



According to the following symmetry relation one can only 
consider J^(l) and J"(2,3) in Eq. i.e. 



°(pqmim2m3i^i 1^21^3 1(1 + P)J\^) 

= 3''(pqTOlTO2TO3Z/lJ/2J^3| 

(1 + P)[J^(1) + J"(2,3)|^), 



(36) 



the matrix elements of the single-nucleon current i.e. J^^(l) 
in the free basis states is independent of the initial and 
final two-body subsystem , and the initial and final spec- 
tators Jacobi momentum vectors are related by the delta 
function (5(q — q' — f Q)- So the matrix element of the 
single-nucleon current takes the following expression with 
the help of Eq. ([16)) and Eq. (l36t and is evaluated as fol- 
lows: 

''(pqmim2m3i/i:/2i^3|(l + P)J^\1') 



3yi .r'" , (q, Q)tf''^r^"' (p,q- -Q) 



'1 ' 13 19 

+ '^m2mi(P - 2^^^)'^n\m„nS-2^ " i^' P " 2^^ ~ 3^^ 
1 

2' 



-KJ''="'\(-P- -q,Q) 

mam' V -t^ o 



7- I^i L'l i-"'? / ^ ^ \ 

X'^m\„nm2(-2P + i"^' " 2^^ " 3^^ 



(37) 



2" 4"' 2" 3 

For two-body current we have: 
°(p'q'm']^m2TO3J^J^'2i^3l J(2, 3)|pqmim2m3i/iz^2i^3)'* 
= S{c{ - q+ iQ)5™,™'/.,.; Jii"r^i:|(p,p'), (38) 
so the two-body current with the help of Eq. p6p is: 
''(pqTOiTO2TO3i/iz/2i^3|(l + P)j"\ip) 

=3 y: fd'p' 



(p, p')^''^''>3 , ( ' _ Iq) 

+ J , , ( p q, v)w / / (p , q 



4 

4q,p)V'„34^'„^(p,q 



4^) 



(39) 



For the meson exchange current we have taken into ac- 
count the momentum transfer of the two nucleons: 



q2 = k2 = + p p', 
q, = k3 kg = iq p + p', 



(40) 



where k2, kj, kg, kg are the initial and final momenta 
of the nucleons 2 and 3. 

Single nucleon current in nonrelativistic limit as well as 
the two-body meson exchange currents are introduced and 
evaluated in the free basis states by suitable expressions 
in the Appendix |21 
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3 Numerical methods and results 

In order to calculate the total cross sections by Eqs. ((ij 
and (HI), we need to obtain the nuclear matrix elements 
using Eqs. (^5)) and In fact, we have to calculate 

the auxiliary state \U) in the free basis stats using Eq. 

To solve this integral equation and also Eqs. ([1]), 
([2]) , (p3)) and ((24)) , we first need to choose the appropriate 
coordinate system. Because of the current properties, We 
have to choose the photon momentum vector, Q, to be 
in z direction, which is also the spin-quantization axis. 
This selection does not impose any restrictions because in 
all of the above equations Q is a constant vector. The p 
and q vectors like any free vectors will be determined with 
three components. Finally, except for the energy of photon 
which is constant, we need six independent variables to 
specify U function in Eq. ([5^ in terms of p and q vectors. 

C/(P, q) = U{P, q, dp, dg, <l)p, (f>g). (41) 

In appendix |B1 The relation between variables of each 
term of Eq. ([5^ and six variables of U function is demon- 
strated. 

For the polar angles which vary from to tt the sign of 
the sine is always positive and we can introduce xe — cos6 
as an independent variable. However, for the azimuthal 
angle, which is in the interval [0, 27r], the sign of the sine 
can not be specified as a function of — cos 4> because 
in each intervals of < (/) < tt or tt < </> < 27r we have 
different signs for sine. So it is difficult to specify the value 
of any function in terms of — cos (f). 

In order to save computing time and computational 
memory we need to define two functions for U in each 
interval. 

U(6)-i^'^'^^ ifO<0<^ 

Ui and U2 are defined as: 

C/i(</>) = F(cos0, - cos20), (43) 
f/2(</>) = F{cos(j), -y/l - COS^(j)), (44) 

We write Eq. (1321) in terms of the independent vari- 
ables of the U function in an appropriate expression for 
numerical iteration. We consider the U operator in a gen- 
eral form: 

\U) = \Uo) + K\U) = \Uo) + K\Uo) + K^\Uo) + ...,(45) 

where \Uo) is the current term and K is an integral kernel, 
then we apply the kernel K to generate the finite Neumann 
series up to (i — l)th order in K . This Neumann series can 

be summed up using Pade approximation to get Up^]^^ 

[35] . We get Upl^^^ from one more iteration. The definition 

of the distance between Up~^^ and f/p^^g is: 



where, the summation runs over all six-dimensional grid 
points. We continue the iteration to reach Z\" < e. e is a 
small number determined by the desired accuracy. In our 
work n = 10 and e = 10"^. 

According to Eqs. (15^ and (flB)) we need the two-body 
t-matrix and triton wave function in the free basis states. 
The extraction of singularity from two-body t-matrix in 
deuteron binding energy (appendix [D| shows that we also 
need the deuteron wave function in the free basis states. 
For calculation of the two-body t-matrix, deuteron and 
triton wave functions we follow the 3D approach intro- 
duced in appendix [C] and we recalculate them in the free 
basis states using AV18 potential [l^. In the calculation 
of the Eq. (15^ we interpolate the calculated data of the 
t-matrix, deuteron and triton wave functions using Cubic- 
Hermit spline method |24) . In the numerical treatment the 
momenta and angles variables should be transformed to 
certain discrete values. For the AV18 potential we use the 
Gaussian quadrature grid points with the hyperbolic map- 
ping for the lower momentum and linear mapping for the 
higher momentum. The numerical Fourier-Bessel transfor- 
mation of the potential encounters difficulties in handling 
at the very high momentum, so it is necessary to use a cut 
off in the integration interval at 150/m~^. 

The meson exchange currents (MECs) were restricted 
to TT- like and p- like exchanges. These MECs are derived 
from AV18 Based on the Riska's recipe [25j . 

We compare the 3D Faddeev calculation, including the 
explicit MECs and the PW representation of Faddeev 
calculation, with Siegert theorem [55]. We have shown 
that the 3D approach with the continuous angle vari- 
ables instead of the discrete angular momentum quantum 
numbers in evaluation of the nuclear matrix elements for 
7 — >■ Nd and 7 — >■ 3A^ leads to less complicated ex- 
pressions but with higher dimensionality of integral equa- 
tions in comparison with the PW representation. 

Calculation of the one-body current and the two-body 
current for a^'^ and are compared together, with 
experimental data in Figs. ([1]) and ([2]). The comparison 
re-confirms the enhancement of the three-nucleon photo- 
disintegration cross section in the peak region and at the 
higher energies. The addition of the three-nucleon force 
can significantly lower the peak in a^'^ and crj*^ and as 
a result gives better agreement with the available data. 
We have displayed in Figs. ^ and (jU that the Siegert 
and MEC predictions are too close in the lower energies 
of the photodisintegration of "^H to Nd and 37V. However, 
the 3D calculation uses the single nucleon current with 
explicit use of tt- and p- like mesons. We found discrep- 
ancies between the two predictions in the higher energies 
of the photodisintegration. We expect that in the higher 
energies of the photon, the use of meson exchange in 3D 
approach produces more sensible results. 

In comparison with the experimental data (Ref |27j). 
Fig. ([3]), the 3D total cross section of — >■ Nd shows 
less agreement at low energies (E^ < 20MeV). The over- 
estimation shows the need for a three-nucleon force (3NF) 
effect. The results of Golak et. al |2H] indicate the improve- 
ment by adding 3NF. The contribution of 3NF in the 3D 
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calculations can be implemented by modifying the \U) as 
an auxiliary state with the term appropriate for the 3NF 
adding to the two-body forces. 



The data of (Ref 



at = 20 - 30MeV due to 
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the insufficient precision can not be compared with the 
theoretical calculations and no concrete conclusions can 
be reached. 

The Skopic et. al data and Kosiake et al data [291130] /\ rnatrix elements of current 

although nearly agree with the 3D calculations in {E~^ = 
15 — 28MeV), however more experimental data is needed 
to reach better conclusion. 

In the case of j^H — >■ 3iV the overestimation of the 
calculated total cross section in the low and medium en- 
ergies in Fig. (4) is also predicted to be due to the ab- 
sence of 3NF. At higher energies more experimental data 
is needed to overcome the discrepancies with the present 
day theories. 



Single nucleon current in nonrelativistic limit consists of 
convection and spin current terms: 



J(l) = GsiQ) 



ki + k; 

2mN 



2m 



N 



GAf(Q)<TX (k'i-ki).(47) 



4 Summary and outlook 



In this paper we have formulated the Faddeev integral 
equations for calculating the two- and three-body pho- 
todisintegration cross sections of the triton in a 3D ap- 
proach. To this aim we have used the free basis states 
which contain Jacobi momentum vectors as well as indi- 
vidual spin and isospin of the nucleons. So we avoid to 
decompose the angular dependent in terms of the angu- 
lar momentum quantum numbers, traditionally used to 
solve these kind of equations, i.e. partial wave approach. 
The final integral equations are less complicated than the 
similar partial wave integral equations and are unique in 
number of the equations in all energies. We have also ex- 
plained how to overcome the moving singularity in Eq. 

by the separation of the radial and angle parts of the 
Dirac delta function. 

Using Eqs. (P) and ^ we have calculated the total 
cross section for 3N and Nd photodisintegration of the 
triton. Benchmarks for the three-nucleon total photodis- 
integration cross sections are presented in Figs. ([3]) and 
©• 

Although the classical approximation of photodisinte- 
gration cross section (predicted by Golak et. al 28 ) and 
the 3D calculation with MFCs are nearly in agreement for 
both 7 Nd and 7 — > 3iV at low energies, the sig- 

nificance of the 3D based calculation with MFCs can be 
tested further with the inclusion of the 3NF calculation. 

Adding the three-body current as well as the three- 
body forces in our calculations are the other major future 
works to be done. We have calculated the two-body t- 
matrices using chiral potential in the 3D approach [17 . 
Therefore the calculation of three-body photodisintegra- 
tion by this potential using the chiral currents is another 
area for consideration. The similar calculation for the ra- 
diative capture is also under preparation. 



Where Ge and Cm are electric and magnetic form factors 
of the nucleon, respectively, ki and k.^ are the initial and 
final momentum of nucleon 1. 

We used the matrix element of the one-body currents 
in the tensor component representation. Considering this 
equality: 



ki -I- k'l = 2q - Q + -(ki +k2 + ka) 
= 2q-Q + ^K, 



(48) 



we can rewrite the single nucleon current in the represen- 
tation of the tensor components as follows: 



jconvec 



TUN 



{Gl{Q)np + Gi{Q)n"), (49) 



J: 



Spin 



V2Q 



2mN 
where 

s+ = \+){+\ 

and finally: 



S±{Glj{QW + Gl,{Q)nn, (50) 



7T" = |n)(n|, 

S- = \-){-l 



(51) 
(52) 



j+'Tnq,Q) 



for V ^ v' 
for m' = TO + 1 



-^^e^'^^GUQ) for to' = TO '(^3) 
-^G-M for to' = TO- 1 



V2 



for V ^ v' 



^Gl,{Q) forTO' = TO + l 



^2)1 



^e-''?'<'G^(Q) for to' 



.(54) 







for to' = TO — 1 



The following tt— and p— meson exchange currents are 
introduced in 9 as follows: 

3^2. = * - GliQ)] (t2 X X3)3 



, (55) 



(T3(T2 ■ q2^'7r(g2) - Cr-2CTs ■ q3W^((?3) 
-1 -yf W^('?3)-l'7r('Z2) )cr2 •q2'T3 ' 
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and 'Vpils) - 



92-93 

Wp(93)cr2 X (0-3 X qg) - i;p(g2)cr3 x ((T2 x Qa) ) ( ^ 2™"!?™/ (g3)i:>+„,„^ (93) ) ) x - q3)±i 

^^p(93) - 'fp(92) 



9 ? 

92-93 



[a, X q,) . K X q3)(q, - q,) -<Z3,±i (^2 m"i^™,™. fe)^+.„j92) 

-0-3 • (q2 X q3)(cr2 X q^) 



-0-2 • (qa X q3)(cr3 x qg) 



(56) 



2 ^ m"D„,'^„,,, (92)£'+"„3 (92) 



The functions w,r((?), Wp(9) and Wp((7) can be extracted ±\/2(q2 ■ q3)'?mf-,,m2±i ( 2 m"£),„f^„t" (g2)£'m"m.,(92) 
from the phenomenological AV18 two-nucleon interaction ™" 

The matrix elements of the two-body current are also +92, ±1 ^2 ^ m D^>_^^„ {Q3)D^i",n2 (93) 
written in the representation of the tensor component. For " 
the pion-exchange we can write: / o \ ^ n / ~ \ n+ f - \ 

(p'q'm'j^m'jmgi/^z^j'^gl | pqm 1 7712^32^1 1^2 ^^3) ™" 

<5™,™;<5.,.;<5(q' - q + ^Q) f G^^(g) - GS(Q)) ^^^^^ ' q3)^™;„™3±i (2 E ^""An^™- (92)I?+.„, (93)) 

(^m^.rns + l ^mj^ ,m3 — 1 )^m2 ,m2 
92x — 93x 



MQs) - vAq2)) -S^'^^^,{6^'^,n3 + l + Sr,^',,m3-l)) (92, - 93,) 



3 

X I (t2\/25™;,,™3±J E m''i:)„^„.(92)£',+ »„,(92) j«^(92) 

±2V2dm'^^m2±l i X! ^"Dm'^m" {h)D^»m3 (93) j ^tt (93) 

(q2 - q3)±i 

92-93 

^^f Xl™"^'"i™"(^2)£'+",„,(92)) +(<5m;,.m2 + l'5™;,.m3-l - S„r'^.r,i2~lS'm[,,r,i3 + l) 

^ 92z - 932 

[5]m"i?„,„„(g3)i?+.„3(93))|. (57) 

In the above equations q2 and q3 are introduced in Eq. 
For the matrix elements of p— exchange current in the (|40l) and Dm'miq) is rotation matrix element for j ~ \^ 
tensor components form we can write: a — 4>q, P = 9q and 7 = 0: 

{]/ c{ 'in\m!2m'^v'^v'2i^'-^\J± |pq™i™2™3'^ii^2^'3) 

A A ' ^^c^^(^P(n^ in~\ Dl,^{aM) = {j^'\R{ah)\3m) (59) 

= '5™,™;<5.,.;<5(q - q + -Q) [G^^iQ) G ^{Q) j ^ ^ 

^ { ^^9i -^91^^^ ('^p(g3) - ^^^(92)^ B U function equation in details 



(58) 



— (=F2A/2(^TO^,,n3±i I E m"D„i'^m„ {q3)D^„^^ (^3) j Wp(q3) The U function Eq. consists of single nucleon current 

V m" ' (SNC), two-nucleon current [TBC) and a complicated 

('^m2,in2+i'^ni3,m3— 1 + '^m2''"2— i'^"i3."i3+i part denoting by(/). 
)(2g3,±it'p(93) - 2q2,±iVp{q2)) 

±2 V2J.,..,, ( X: m"D,^,^„ iq2)Di„^^ (^2)) .p(92) + ' ^^^'^ (60) 
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For the SNC part we can incorporate the one-nucleon sin(j)2p 
current as follows: 



SNC 



3E (9.^.>^^jc\i?™3(^''5''^p'^«'^^.'^^.) (^1) 



+ (^2, X2q, X24,^ )i'2lmsm, (P2, 92, X2p, 4,,, X24,^ , aJj^, ) 



+'^±,m3mi(«3,a;3<„a;3^,) 

, I/( L^l 1/2 / / / \ 

X V'^^mim^ (^'3, 93, ^Sp, X^g, X34,, , X^^^ ) 



X2<t,g = COS (j)2q = 



pJl - xlx^^ - ^qJl - xlx^ 



(61) 



g2\/l - x"^ 



(72) 



2q 



In term of the p, q, Xp,Xq, x^^ , x^^ we can write all the 
variable as follows: 

q' = \<l-lQ\ = {q' + lQ^-lqQxq)K (62) 



13 19 3 1 
P2 = \- ^P- ^q.\ = (^/ + Y^q^ + -pqcos^)^, (63) 



92 = |P - ^q| = (/ + + pqcosj) = , (64) 



cos 7 = XpXg + Y^l - x'^^1 - X; 

X<f,pX^^ 



i-xi ^ 1- X 



(65) 



sin ^2,= (p^Ji - xl^ 



/q2^l-x\ 



(73) 



P3 



93 



Jq| = (J/ + ^9^ - ^P9cos7)i, (74) 



P-^q| = (/ + j9^+P9cos7)^ (75) 



X^q = 



93 



(76) 



X2q 



pXp 



2 9^9 



92 



(66) 



X^p 



1 3 

^pa^p + -^qXq 

P2 



(77) 



a:2p 



P2 



92 =|P-2q-3QI 

1 2 4 
4^+9 
4 2 a 



(P^ + 79^ + -P9C0S7 



2^2g — 



pXp 2 9a'(/ 3 Q 



9^ 



(67) 



(68) 



(69) 



1 2^, 
93 = |-P - 3QI 

= (P^ + ^9^ + -|-pgcos7 

4 2 1 

+-^pQxp + -^qQxq)^, 



X3<l>p = COS ^3p 

= (-5.\/^=^*. 
/(k\/i 



»V^-*.) 



' 1 ^3p / ' 



(78) 



(79) 



a;20p = COS ( 



'2P\J^ XpX<Pp 49yl 2;2a;0^ 



P2^1 - X 



(70) 



sin (j)3p 



-p^l - x'p^l - x^2 + ^9^1 - -'4^/1 - xj 



/P2^ 



^-4p- 



(80) 
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= cos 03^ 



/( q2Jl-xl„ ), 



(81) 



g2 \/l - 



3q 



(82) 



For two-body current, i.e. the second term in Eq. (|60p 
one can write as follows: 

TBC = 3 Yl r d^P'p'^ f ^^M= 



Ftbc{x4,,,JI 



+Ftbc{~x^^,,-^II - x-^J \, 
Where 



(83) 



Ftbc{x4,^,, ^l-xl^^) 

X Ci m?m^ {p',q',Xp,,Xq,,X^^,,X^^) 
X CTm?m^ (^'' ,q',Xp',Xq,,X^^,,X^J 



(84) 



sin (/i2b = (^p^ 1 - a;2 y' 1 - a;|^ - 1 - 1 



92b 



+4? +P 



/2 



926 



= l2Q + p-p'l 

= i^Q^ + P^ + p'^ + pQxp 
-p'Qx'p - pp' cos /3)i, 



cos 13 — Xpx'p + y/l - .x^ y^l - x'^ 



X(ji^,x^^ + \/l xf. 



92 



(85) 



(86) 



(87) 



3 1 
-f-gp' cos /?')*, 



cos/3' = Xp'Xq + y'l - 

X(j,X(j) , 



1 - a;i,/l - a; 



^2fc 



-^9\/l -^g^;^, -P' \J^-x'^x^^ 



In term of p,q,Xp,Xq,x^^,x^^ we write the j's argument sin02h 
as follows: 



qJl-X^Jl~xl-p\l-x'\l 



926 



/ 9i^/l-^n> 



l-Q -p+-q p' 



1 



-^Qgzg -f pp' cosf3 — ^pqcos7 

3 1 
-qp' COS /3')^i 
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'''2b 



smc 



-q\ll- xlx4,^-^ x'^x^^, 



(97) 



qJl-xlJl-xi+c('Jl-xi,„Jl-x" 



/(l2q + q"IVi-4 



3q^ 3g2 

4 4' 



(106) 



(107) 



-aJl^xlJl-xl^-p'Jl-x';^ 



l-x2 

0p 



(98) 



- 2 ^ 

lh" + q| ^ 



(108) 



|q+^q"l = (j9"' + 9" + 99'V')*, (109) 



g' = |q - ^Ql - (9' + - \qQxq)"^, 



Eg' 



([Xq ^ 

7^' 



(99) 



(100) 



XqXqii + ^1 ^qy^ ~ ^q" 
Xcj)"X^^ 



l-xl„Jl-xl 



The third part of Eq. (1601) in term of the integration vari- 
ables and the six variables of the U function can be written 



as: 



For delta function we have: 

5{x" -xo)^6[f{xq„)], 

where: 



(110) 



(111) 



1 



qE-Ed 



+p 

\i-p'\ 



q"q" 



3q 



dx„ 



E 



dp 





dx'l 



E — 3«i _ El 

4m m 



f(Xqn) = aXqn + \/ 1 - O? ^ I ~ x'^^„C - Xq, (112) 

in the above equation, a — Xq, c — x^ ,, and xq = 



-1 ^n-x^„ 



qq" 



X < F{x^n , y'l - xl„ ) + F(-X0" , - 1 - ) 
x5{x"-xo), (101) 
Where 



For handling the Eq. (|llll) we have to use delta func- 
tion properties to change the argument from x" to Xqii . 
So we need zero points of f{x), i.e. Xi: 



(113) 



c2(a2 — 1) — ' 



f^(x^.- , ^1 - x2 „ ) = (p, p' , , Xp- , , x^^, ) Where 

xC/(p",g",Xp",x,",x^^„,,, J, (102) A = ^{a^ - l)[c^{a' - 1) + - a^] 



(114) 



iqa:, -I- q"xqn 

lh + q"l 



(103) 



C Two-body t-matrix, Deuteron and Triton 
wave functions in free basis states 



-, , ^ Two-body calculations in the 3D approach are performed 

|-q-^q"| = {-q^ + q""^ + qq"x")^ , (104) in Refs. We briefly introduce the two-body t- 

^ ^ matrix in the momentum-helicity basis states as: 

X., = cos <i>^, ''%WSAT\t\^vSATY- = t\^A^{p\ p), (115) 



~ XqXcjjg + q" ^ Xq"x'd, Where the momentum-helicity basis states is defined by: 

— ^= , (105) 

|iq + q"L/l-4 \ppSATr'^ = {l-fj^{-lf+^)\ppSA)^\T), (116) 
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where S is spin, A is its projection to p and T is the isospin 
of the two-body system. If we choose the direction of the 
vector p to be in the z-direction we can write: 



(117) 



and the Lippmann-Schwinger equation in the momentum- 
helicity basis states for any energy is written as: 

tl,^/{p',p,Xp,) = VX,r{p\p,Xp,) 
.1 



4i: 



112 

1 q 



xG^{Ep)t-JJ{p'\p,x"), 



1 II ttST.A/ I II ii\ 

dx Vjy,^,,' (p ,p ,Xp',x ) 



(118) 



Where 



IT ST, A, I II iix 



b"e~^^(^'-^")VXd^{p',p"). (119) 



The relation between t-matrix in the momentum-helicity 
basis states and those in the free basis states is: 

" (p' m'j I i I pm2 ma i/2 ) " 

1 „ 



^(1^2+1^3), ( 

AST 



'22 
,1 1 



'22' 



'22 



X 



'22 

j-ttST (I 



AA' 



tA'Y(p',P), 



(120) 



rotation matrix element, (if,^(a;p) is introduced in (I59p . 
The deuteron wave function in the free basis states is also 
related to the momentum-helicity basis states: 



(V'^-'^''|pm2m3i/2J^3)'' 
_ 1 
^ 2 



+4>r(p)d\,,{xp) - 4>"'{-v)d\,^i{xp) 
xC(iil,m2m3ylo)C(iio,i^2J^30) 

piAo4>qjl I \ 



(121) 



Where Md is spin projection of the deuteron, we have 
summed over A = —1,0,1, and we have used this fact 
that the terms with A = —1 and A = 1 are the same. Also 
we define: 



la / 



{qqlA;0\i//'^) = ^^'^ {q, x,)e^^'-^^ , (122) 



the deuteron wave function in the momentum-helicity ba- 
sis states is calculated using eigenvalue equation: 



{^-Ed)^^'-{q,X,) + ^ 



dqq^ 



1 C°° 

dx'v\f-{q,q',x,,x')cP^'-iq',x') + - / dqq' 
' dx'v\l'-''''{q,q',Xg,x')cf>^'-{q',x') = 0. (123) 



The triton wave function is calculated by applying the 
formalism of Ref.|16|. In this formalism the triton wave 
function has been evaluated in the basis states: 

{pqa\W) = (pq(s23, \)SMsit23, \)TMt\^), (124) 

In the above equation S23 is total spin of the subsystem 
5 and Ms are the total spin of the three particles and its 
projection along the z axis respectively. The same expla- 
nation is used for the isospin. 

The triton wave function in Eq. (jl24p is reproduced 
using AV18 potential and is related to the one in the free 
basis states as follows: 

{pc\mim2m^vii'2V7i\^) = y^{mim2m3iyii'2i^3\a){pqa\^) 

a 

= Y.9^o.{p(ia\^). (125) 

a 

\a) and I7) are spin and isospin states that are defined by 
Eqs. (jl26p and (jl27p and the Clebsch-Gordan coefficients, 
g^a, are introduced in Ref. [16) . 



|a) = I(S23, l)SMs{t23, \)TMt), (126) 



I7) = \mim2m3i^ii^2V3) ■ (127) 



D Extraction of Deuteron singularity 

To extract singularity of the two-body t-matrix in deuteron 
binding energy we need to evaluate the following term in 
the momentum-helicity basis states. 



lim {E-Ed)t\Z 

E^Ed 

= ^''{p'p'lA';0\V\W^^''){^l^''''\V\pplA;Oy'', (128) 

two terms on the right hand side of the above equation 
evaluate separately. 

i'^(p9lyl';0|T/|iff'') 

= iE / rfp"'"(p9i^';0|v^|pV"i^";0)i'^ 

A" 

X i''(pVi^";0|tf'*^'') 

= iE / dp"V}}Up',P")1>A'fip'% (129) 
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and 



[ dp"{^^'"\V\p"p"lAOy'' 

A" 



X -{p";p"lAO\V\pplAOy 



la /„" 
1 

A' 



4E 



9. J.Golak et al. , Phys.Rept. 415, 89 (2005). 

10. R. A. Rice, Y. E. Kim, Few-Body Syst. 14, 127 (1993). 

11. I. Fachruddin, Ch. Elster, W. Glockle, Phys. Rev. C 62, 
044002 (2000). 

12. I. Fachruddin, Ch. Elster, W. Glockle, Phys. Rev. C 63, 
054003 (2001). 

13. M. R. Hadizadeh and S. Bayegan, Eur. Phys. J. A 36 
201(2008). 

dp"V\}P^{p" ,p)<P*J^!'' {p"). (130) 14. S. Bayegan, M. R. Hadizadeh, and W. Glockle, Prog. 

Theor. Phys. 120 887 (2008). 
15. S. Bayegan, M. Harzchi and M. R. Hadizadeh, Nucl. Phys. 
A 814 21(2008). 



By multiplying these two terms and considering the 



derive the following equation using Eq. (jll9l) 
lini iE~Ed)t\Zip',P,x') 

E^Ed 



p"^dp"dx" 



110, A/d 



{p',x\p",x")K'{p",x") 



+2v]l,l^^'^{p\x\p",x")<Pr'{p",x") 



momentum vector p to be in z direction,we can extract the 16. S. Bayegan, M. R. Hadizadeh, and M. Harzchi, Phys. Rev. 
azimuthal angle parts of the potential and the deuteron C 77, 064005 (2008). 

wave function by applying Eqs. PTT)) and We then 17. S. Bayegan, M. A. Shalchi, M. R. Hadizadeh, Phys. Rev. 
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p"'dp"dx" [VoTip'\p, x")1>l''^ ip", x") 



-f2y/7(p",p,x'>-"''(p",:r") 
1 



- I p"^dp"dx"^ 



(^v)lP,\p\x\p\x")1>{^{p'\x") 
+2v)lP,^{p', x',p", x")^>tip", x")) e-'*^^(*'-*" 
p"^dp"dx" (v„'rip'\p, x")1>l\p\x") 
+2V^\\p\p,x")<!>l\p",x")^ 



(131) 



In above equation we have used this equality: 



d^"e^(A-A/.)(f"-<i.') ^ 2T,5AMd 



(132) 
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Fig. 1. Comparison of the 3D calculations for the total cross section of the two-body (Nd) photodisintegration of Triton using 
the single nucleon current (dashed line) and the two-body current (solid line). The experimental values are taken from Ref. [27) 
as EXP. DATAl, [29] as EXP. DATA2 and [30] as EXP. DATA3. 




Fig. 2. Comparison of the 3D calculations for the total cross section of the three-body (3N) photodisintegration of Triton using 
the single nucleon current (dashed line) the two-body current (solid line). The experimental values are taken from Ref. [27) . 
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Fig. 3. Total cross section for Nd photodisintegration of Triton. The solid line is the 3D calculation with AV18 potential. The 
dashed line represents the partial wave calculation using AV18 potential and dot-dashed line represents the partial wave result 
using AV18 potential and 3NF (Urbana IX). The experimental values are the same as Fig. ([l}. 
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Fig. 4. Total cross section for 3N photodisintegration of Triton. The solid line is the 3D calculation with AV18 potential. The 
dashed line represents the partial wave calculation using the AVIS potential and dot-dashed line represents the partial wave 
result using AVIS potential and 3NF (Urbana IX). The experimental values are the same as Fig. 



